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I present in this paper the neutron spin physics program in Hall A at Jefferson Lab- 
oratory using a polarized 3 He target. The program encompasses several completed ex- 
CNj ■ periments, in which, valuable spin observables (spin dependent structure functions) were 
measured in order to learn about how the nucleon spin arises from the behavior of the 
constituents. These experiments also offer a ground for testing our understanding of 
the strong regime of quantum chromodynamics (QCD) the theory of strong interactions 
through the determination of moments of these structure functions. 

x ■ 

1. INTRODUCTION 

<d : 

Technical advances for producing polarized beams and polarized targets in the early 80 's 
triggered a novel experimental effort at SLAC which provided for the first measurement 
of the proton spin asymmetry A±(x,Q 2 ) PQ. The limited x range of this measurement 
led to an interpretation of the results consistent with the predictions of a naive non- 
relativistic constituent quark model (CQM). However after EMC at CERN [2] extended 
the range of the SLAC measurement to significantly lower values of x, a different and 
surprising picture of the spin structure of the nucleon emerged. By combining the mea- 
sured Ti = Jq gi(x 1 Q 2 )dx with the neutron and hyperons (3 decays within the CQM 
it was found that only a small fraction of the total nucleon spin was accounted for by 
quarks. Since then, enormous theoretical and experimental progress has been achieved 
in understanding the nucleon spin structure |3|4j . The original findings of EMC were con- 
firmed and measurements on the neutron allowed the test of the Bjorken sum rule [5] to 
< 10% level. This accurate test was only possible after the sum rule, originally derived 
at Q 2 — > oo using current algebra, was re-derived using the technique of operator prod- 
uct expansion (OPE) within QCD |6|7|8|9j . This step was essential in order to calculate 
the corrections necessary to evolve this sum rule to Q 2 values accessible experimentally. 
In these studies most of the measurements on the nucleon spin structure functions were 
performed at large Q 2 above 1 GeV 2 . Nevertheless to take full advantage of the OPE and 
test its applicability in the strong regime of QCD, it was understood that measurements 
of the nucleon spin structure below Q 2 = 1 GeV 2 and in the resonance region would be 
important for determining the size of the higher twists corrections. In this connection it 
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was realized that the Bjorken sum rule was only a special case of a more general sum rule 
known as the "extended" Gerasimov-Drell-Hearn (GDH) sum rule which coincides with 
the GDH sum rule at Q 2 = and with the Bjorken sum rule at Q 2 — > oo. 

Furthermore, higher moments of the spin structure functions are connected with specific 
twist matrix elements (observables) which can be evaluated using Lattice QCD |lU|llj . 
The advantage of measuring higher moments of the spin structure functions is twofold, 1) 
the kinematical region which gives most of the contribution to these moments is experi- 
mentally accessible 2) calculations of specific matrix elements related to these moments 
through OPE are possible using Lattice QCD. 

I will present here results of two experiments performed at Jefferson Lab in Hall A that 
address two areas of the above mentionned concerning the spin structure of the neutron. 
First we present results of a precision measurement (E99-117) of the neutron asymmetry 
A™ in the valence region where the existing world data could not confirm the expected 
behavior from either the CQM or pQCD due to a poor statistical precision. Second we 
present results of a measurement of both g\ and gi (E94-010) through the resonance region 
at a Q 2 ranging from 1 GeV/c to 0.1 GeV/c and the neutron T™ and are evaluated in 
the same range of Q 2 connecting the perturbative regime to the strong regime of QCD. 



2. EXPERIMENTAL METHOD AND SETUP 

The experiments described in this presentation share the same experimental setup |12j . 
They were carried out at Jefferson Lab in Hall A using a highly polarized electron beam 
(70-80%) with an average current up to 15 n A and a high pressure polarized (on average 
between 30% and 40% in-beam) 3 He target with the highest polarized luminosity in the 
world. The target is based on the spin exchange principle where rubidium atoms are 
polarized by optical pumping and their polarization is transferred to the 3 He atoms via 
spin exchange collisions ^J). The target was polarized either parallel or perpendicular to 
the direction of the electron beam and its polarization measured by three independent 
methods: NMR through the Adiabatic Fast Passage technique, Electron Paramagnetic 
Resonance technique (EPR) and elastic scattering off 3 He. In each case the results of the 
three methods agreed within systematic errors. 

The scattered electrons were detected using two HRS spectrometers which sat symmet- 
rically on each side of the electron beam line to double the count rate. The spectrometers 
were equipped with standard detector packages which consists of two vertical drift cham- 
bers for momentum and scattering angle determination, a CO2 gas Cherenkov counter 
and a double-layered lead-glass shower counter for particle identification The elec- 
tron identification efficiency was 99% and the 7r~ rejection factor was found to be better 
than 10 4 for both spectrometers which was sufficient for all experiments. The momentum 
of each spectrometer was stepped through the relevant excitation spectrum at each inci- 
dent energy of a given experiment. In many cases the comparison of the data between 
the two identical spectrometers at symmetric angles allowed us to minimize many of the 
systematic uncertainties. 
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3. SPIN AND FLAVOR DECOMPOSITION IN THE VALENCE QUARK 
REGION 

Among the deep inelastic lepton scattering (DIS) spin observables important for under- 
standing the structure of the nucleon at large x, the virtual photon-nucleon asymmetry 
A™ and spin structure function are the most poorly known. This shortcoming is due 
to the small scattering cross sections in the large x region combined with the lack of high 
polarized luminosity facilities. In the past experimental effort has gone into measuring 
A™ and g™ at low x, since it provides for most of the strength to r 1; resulting in a poor 
statistical precision in the world data for x > 0.4. This region, however, is clean and 
unambiguous since it is not polluted by sea quarks and gluons and thus offers a unique 
opportunity to test predictions that are usually difficult if not impossible at low x. A very 
clean contribution of the "valence quarks" can be expected when x > 0.5. At present at 
large x the world data are still consistent with the most naive of CQM 's which uses a static 
SU(6) symmetric wave function to describe the ground state of the nucleon. However, we 
know from the ratio of neutron to proton unpolarized structure functions measurements 
that this symmetry is broken in nature and thus the prediction is unrealistic. 

The set of predictions of A™ in the valence quark region falls into two categories, those 
of CQM's which break SU(6) symmetry in the ground state wave function by hyperfine 
interaction |15|16j which responsible for the N-A mass splitting, and those of perturbative 
QCD with a hadron helicity conservation (HHC) constraint |17|18j as x — > 1 which break 
SU(6) symmetry dynamically. While both classes of models predict that A" — > 1 when 
x 1 their x dependence is noticeably different, allowing for a softer variation of A™ 
towards unity in the case of the CQM's. In the case where the HHC constraint is not 
imposed in pQCD calculations the results are rather a fit to the world data (201 ■ m 
this case the approach is to perform a pQCD analysis of the proton and neutron world 
data of the polarized structure functions at next-to-leading-order (NLO) and produce a 
global fit with a minimal number of parameters . This endeavor has been successful in 
the calculations by Bourelly, Soffer and Bucella ,21, using a statistical physical picture 
of the nucleon. 

The difference between the approaches described above is more pronounced when the 
constituents flavor decomposition is performed. In the case where we consider a proton, 
we have Au(x)/u(x) — > 1, Ad(x)/d(x) — > 1, for the case of pQCD with HHC models, and 
Au/u -> 1, Ad/d -> -2/3 for the case of CQM's. We notice that in the pQCD with HHC 
models Ad/d changes sign from negative at low x to positive at large x. This behavior 
is clearly different in the CQM model based on single flavor dominance where Ad/d is 
negative at all x values. 

Furthermore, we mention exceptions to the models discussed above. A bosonized chiral 
quark model where the nucleon appears as a soliton in the meson fields was used to 
calculate A\ in [22] • In !2H] local quark-hadron duality is used to link the structure 
functions g\ and F% at large x with the Q 2 dependence of the nucleon electromagnetic 
form factors. A prediction was possible however we do not know how low in x will this 
prediction works. Finally structure functions have been calculated in the bag model and 
the effect of the pion cloud on these functions was evaluated in |24j . 

At Jefferson Lab we are in a unique position to take advantage of the unprecedented 
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polarized luminosity. The experiment was carried the highest available incident energy 
of 5.7 GeV. Both the An and Aj_ asymmetries for 3 He were extracted after correcting the 
raw asymmetries by the beam and target polarizations and the dilution factor at three 
kinematics: x = 0.331, 0.474 and 0.61, and Q 2 = 2.738, 3.567 and 4.887 GeV 2 , with W 2 = 
6.426, 4.846 and 4.023 GeV 2 respectively. False asymmetries were checked by measuring 
asymmetries of the polarized beam scattering off an unpolarized 12 C target and found 
to be negligible. Radiative corrections were applied to obtain the Born asymmetries and 
then, the asymmetry A\ and the ratio of polarized to unpolarized structure functions 
gi/Fi of 3 He were extracted using 
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where D, d, r\ and £ and D' are kinematical variables. Finally a 3 He model which includes 
the S, S' and D states and pre-existing A components in the 3 He ground state wave 
function was used to extract A™ and gi/iqof the neutron from that of 3 He [26 . 

In Figure 1 (right panel) we show preliminary results of A™. The first data point at 
x = 0.33 is in good agreement with previous measurements. The data points show a 
clear change of sign of A™ as x increases and are compared with theoretical predictions. 
The total error in each point is dominated by the statistical error. The solid line is a 
pQCD HHC based on LSS(BBS) parameterization of gi/F™ [TP] , the long-dashed line 
is a prediction of gi/F™ from LSS 2001 parametrization at Q 2 = 5 GeV 2 without HHC 
constraints 20J. The shaded area is a range of predictions of A™ from the constituent 
quark model ^H] while the dot-dashed line is a calculation of the statistical model at Q 2 
= 4 GeV 2 by Bourrely et. al. |2T], and the dotted line is the local duality prediction by 
Melnitchouk (2H|. Finally, the short dashed line is the chiral soliton model prediction at 
Q 2 = 3 GeV 2 by Weigel, Gamberg and Reinhardt [22J. Data from Hermes and SLAC are 
original values without being re-analyzed for the A contribution of the nuclear corrections. 

Assuming that the strange quark distributions are negligible in the region 0.3 < x < 1 
we used the quark parton model interpretation of gi and F\ to perform a flavor decom- 
position of the spin dependent quark distributions. 



Au _ 4 g{ f d 
~u ~ 15Ff V + u 



1 gl A A d 
— — 1 + 4- 
15 Ff \ u 



d 15 F? V d 



Here u and d represent the sum of quark and antiquark distributions. Using equation (2) 
and the d/u ratio extracted from the proton and deuteron structure functions data 
we present in Fig. 1 (left panel) results of the down quark distributions obtained in 
E99-117 (filled squares) along with preliminary results of the HERMES semi-inclusive 
measurements (open circles) j2Zj. The solid line is a pQCD fit to the world data using the 
HHC constraint as x — > 1. The dashed line correspond to a non relativistic constituent 
quark prediction. It is clear that up to x = 0.6 the data favor the CQM and is in violation 
with the HHC pQCD based calculations. We point out that this is consistent with the 
interpretation the recent Q 2 dependence of the proton electromagnetic form factors ratio 
measured at Jefferson Lab 
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4. EXTENDED GDH SUM RULE AND TWIST-THREE MATRIX ELE- 
MENT 

In its extended form the GDH sum rule spans all range of momentum transfers and 
reads |30 a : 

4 /,, 2 dv = S 1 (Q ), 3 

where gi(v, Q 2 ) is the spin structure function of the nucleon, Si the virtual Compton scat- 
tering amplitude, v the energy transfer and M the nucleon mass. Starting the integration 
in equation (3) from the pion production threshold this sum rule can be expanded around 
Q 2 = using for example the heavy baryon chiral perturbation theory (HB^PT) |31j . 



r?(Q 2 



Q 2 2 5.54 



8M 2 



2M 2 



Q 4 



(4) 



where K n is the anomalous magnetic moment of the neutron. The slope of the first term 
in the right hand side of equation (4) corresponds to the standard GDH sum rule [37] 
while the second term on the right hand side has been evaluated using HBxPT. Other 
evaluations are available [35 3Qj and will be compared to the data below. In the large 
Q 2 region Eq. (J3J takes another form. If we consider three quarks flavors and three-loop 
result for the twist-two component r™(Q 2 )is given by [32J: 
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where g a = 1.257 and a 8 = 0.579 [o3j are the triplet and octet axial charge, respectively. 

is defined as the renormalization group invariant nucleon matrix element of the singlet 
axial current [33]. At large enough Q 2 the /14 contribution becomes negligible. 

Presently, there is a large set of data on the spin structure function at Q 2 > 1 GeV 2 
used to determine and a s . However, below Q 2 = 1 GeV 2 the experimental situation 
is less than ideal, prompting for the investigation of higher twists as Q 2 decreases below 
1 GeV 2 . This experimental situation has now changed with the completion of experiment 
E94-010 which covers a Q 2 ranging from 0.1 GeV 2 to 1 GeV 2 . In this experiment we have 
measured the helicity dependent electron scattering cross sections for excitations energies 
covering from the quasielastic region up to the deep inelastic region, thus including the 
resonance region. The longitudinal (Aerii) and transverse (AcxjJ cross section differences 
are formed from combining data taken with opposite beam helicity, 
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ry = a- -a--, Affi = a 1 ^ - cr T ^, (6) 

where o-WW and cr^(T=>) are S pj n dependent inclusive differential cross sections with 
the electron spin helicity anti-parallel (parallel) and anti-perpendicular (perpendicular) 
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to the target spin direction. The 3 He spin structure functions gi(x, Q 2 ) and g 2 {x, Q 2 ) were 
extracted from the expression: 



gi(x,Q 2 ) 
92(x,Q 2 ) 



MQ 2 v E 



Aa 2 e E'E + E' 
MQ 2 v 2 1 



4a 2 2E'(E + E' 



Aery + tan(0/2) Aa ± 



E + E' cos 9 . 
11 E< sin 9 



(7) 
(8) 



where i±7 and E' the incident and the scattered electron energy respectively, 9 is the scat- 
tering angle, M is the nucleon mass and a e is the electromagnetic fine coupling constant. 

The integral Fi(Q 2 ) for the neutron is evaluated by first performing the integration 
of the 3 He response at constant Q 2 from the nucleon single pion production threshold 
to the lowest x measured in the experiment. We then subtracted a small contribution 
due to the quasielastic tail and applied nuclear corrections to the integral following the 
procedure described in ref. [38J. We estimated the missing deep inelastic contribution in 
the range 2 GeV < W < 30 GeV according to the parametrization of Thomas and Bianchi 
parametrization [39] . 

In Fig. El (left panel), we show T"(Q 2 ) at five values of Q 2 from 0.1 to about 1 GeV 2 
by step of 0.2 GeV 2 (solid circles) compared to the world data (open symbols). The 
band located around zero represents the size of systematic errors. At low Q 2 the lines 
correspond to calculations by Bernard et. al. [US] (short dashed-dot line) without 
vector mesons and by Bernard et. al. [HB] (dotted line) with vector mesons, respectively. 



The long dashed line is a calculation by Ji et. al. [30131] using heavy baryon x?T. The 
solid black line represents the GDH result. The comparison with our data shows that 
beyond Q 2 =0.1 GeV 2 these calculations do not reproduce the trend of the data because 
the A degrees of freedom might become important but are not included in both 
calculations. At moderate and large Q 2 two calculations, one by Soffer and Terayev [12] 
(short dashed curve), the other by Burkert and Ioffe (13] (long dash-dot line) are plotted. 
Soffer and Terayev assume that the integral over g? = gi+g2 varies smoothly from high Q 2 
where #2 ~ down to Q 2 = 0. Using their prediction for this integral and subtracting the 
contribution from #2 using the Burkhardt Cottingham sum rule 01] gives the short dashed 
curve in Fig. El which agree relatively well with the data. Burkert and Ioffe consider the 
contributions from the resonances using the code AO, and the nonresonant contributions 
using a simple higher-twist type form fitted to the deep-inelastic data. Their model is 
constrained to fit both the GDH and the deep-inelastic limits, and it describes the data 
quite well. 

In the high Q 2 regime, using OPE, the di matrix element defined as 

d 2 (Q 2 ) = [ x 2 [2 9l (x, Q 2 ) + 3g 2 (x, Q 2 )}dx, (9) 
Jo 

is shown to be a measure of the electric and magnetic polarizabilities of the color field [3]. 
At small Q 2 , a region covered by our data, its conventional interpretation in terms of 
higher twist is not obvious. A recent publication [3T] gives a new insight in this regard at 
low Q 2 using rather xPT. Thus, the Q 2 evolution of o? 2 is also a quantity that can shed 
some light on the strong interaction in the nucleon. 
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In Fig. |21 (right panel), the c?2 matrix element is shown at several values of Q 2 where the 
integration in equation (9) excludes the elastic peak. The results of this experiment are the 
solid circles and the grey band represents their corresponding systematic uncertainty. The 
SLAC E155 [46J proton (open circle) and neutron (open square) results are also shown. 
The solid line is the MAID calculation 40] while the dashed line is a Heavy Baryon x?T 
calculation j?T] valid only a very low Q 2 . The Lattice prediction [TT] at Q 2 = 5 GeV 2 for 
the neutron matrix element ( not shown here ) is negative but close to zero. We note 
that all models predict d\ to be negative or zero at large Q 2 . At moderate Q 2 the data of 
E94-010 show a positive d% but decreasing perhaps to zero at high Q 2 . The SLAC data 
also show a positive d% value but with a rather large error bar. More measurements are 
needed to have a complete determination of the transition from low to very high Q 2 of 
this important matrix element. 

5. CONCLUSION 

In summary, we presented results of two experiments at Jefferson Lab Hall A which 
took advantage of the highly polarized beam and high pressure polarized 3 He target 
to investigate the internal spin structure of the neutron in the perturbative and the 
strong regimes of QCD. In E99-117 we have determined the world most precise spin 
and flavor dependent distribution in the valence region. The results are more in line 
with the constituent quark model than with the HHC constrained pQCD prediction. In 
experiment E94-010 we have measured both gi and #2 m the resonance region allowing 
for the first time to study the neutron integrals and the twist three matrix element d"% 
from the pQCD regime to a regime where chiral perturbation theory might give us more 
insight into the structure of the nucleon. With time we expect Lattice QCD to make 
predictions in the intermediate region. 
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